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Abstract 

Maxwell's equations are introduced in their general form, together with a basic 
set of mathematical operations needed to work with them. After simplifying 
and adapting the equations for application to radio frequency problems, we 
derive the most important formulae and characteristic quantities for cavities 
and waveguides. Several practical examples are given to demonstrate the use 
of the derived equations and to explain the importance of the most common 
figures of merit. 



1 Introduction to Maxwell's equations 
1.1 Maxwell's equations 

Maxwell's equations were published in their earliest form in 1861-1862 in a paper entitled "On physi- 
cal lines of force" by the Scottish physicist and mathematician James Clerk Maxwell. They represent a 
uniquely complete set of equations that covers all areas of electrostatic and magnetostatic problems, as 
well as electrodynamic problems, of which radio frequency (RF) engineering is only a subset. Surpris- 
ingly, they include the effects of relativity even though they were conceived much earlier than Einstein's 
theories. 

In differential form. Maxwell's equations can be written as 



VxH = J + ^, (1) 
dB 

V X E = — ^ , Maxwell's equations (2) 

at 

V-D = ^v, (3) 
V-B = 0, (4) 

where the field components and constants are defined as follows: 

E, electric field (V/m) 

D = eogrE , dielectric displacement (A s/m^) (5) 

B , magnetic induction, magnetic flux density (T) 

H = , magnetic field strength or field intensity (A/m) (6) 

J = kE, electric current density (A/m^) (7) 

d , 

— D , displacement current (A/nr) 
at 

eo = 8 . 854 • 1 0" , electric field constant (F/m) (8) 

gf , relative dielectric constant 

/io = 4;r • 10~^ , magnetic field constant (H/m) (9) 

, relative permeability constant 

K . electrical conductivity (S/m) 



In the following sections, we shall see that most of the important RF formulae can be derived in a 
few lines from Maxwell's equations. 



1.2 Basic vector analysis and its application to Maxwell's equations 

In order to make efficient use of Maxwell's equations, some basic vector analysis is needed, which is 
introduced in this section. More detailed introductions can be found in a number of textbooks, such as 
for instance the excellent Feynman Lectures on Physics |[T|. 



Gradient of a potential 

The gradient of a potential is the derivative of the potential function ^{x,y,z) in all directions of a 
particular coordinate system (e.g., x, y, z). The result is a vector that tells us how much the potential 
changes in different directions. Applied to the geographical profile of a mountain landscape, the gradient 
describes the slope of the landscape in all directions. The mathematical sign that is used for the gradient 
of a potential is the 'nabla operator'; applied to a Cartesian coordinate system, one can write 
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gradient of a potential 



(10) 



The expressions for the gradient in cylindrical and spherical coordinate systems are given in Appendices 
EandO 



Divergence of a vector field 

The divergence of a vector field a tells us if the vector field has a source. If the resulting scalar expression 
is zero, we have a 'source-free' vector field, as in the case of the magnetic field. From basic physics, 
we know that there are no magnetic monopoles, which is why magnetic field Unes are always closed. 
In Maxwell's equations Q, this property is included by means of the fact that the divergence of the 
magnetic induction B equals zero. 

In Cartesian coordinates, the divergence of a vector field is defined as 
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divergence of a vector 



(11) 



The expressions for cylindrical and spherical coordinate systems are given in Appendices [B] and [C] 



Curl of a vector field 

When we form the curl of a vector, we are interested in knowing if there are any curls or eddies in the 
field. Let us imagine that we are looking at the flow of water in a cooling pipe. To check for curls, we 
can use a stick around which a ball can rotate freely. We position a Cartesian coordinate system at an 
arbitrary origin and align the stick first with the x axis, and then with the y and z axes. If the ball starts 
rotating in any of these positions, then we know that the curl of the vector field describing the water 
flow is non-zero in the direction of the respective axis. The curl of a vector a is therefore also a vector, 
because its information is direction-specific. Its mathematical form in Cartesian coordinates is defined 
as 
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curl of a vector 



(12) 



The unit vectors u„ have no physical meaning and simply point in the x, y, and z directions. They have 
a constant length of 1. The expressions for cylindrical and spherical coordinate systems can be found in 
Appendices [B] and |C] 



Second derivatives 

In some instances, we have to make use of second derivatives. One of the expressions that is used 
regularly in electrodynamics is the Laplace operator A = V^, which — since the operator itself is scalar — 
can be applied to both scalar fields and vector fields: 

A0 = V • (V0) = V2(/) = 44 + ^ + ^ • Laplace operator (13) 

The expressions for cylindrical and spherical coordinate systems can be found in Appendices [B] and [C| 
We also introduce two interesting identities: 

Vx(V0)=O, (14) 
V-(Vxa)=0. (15) 

Equation ( [141 ) teHs us that if the curl of a vector equals zero, then this vector can be written as the 
gradient of a potential. This feature can save us a lot of writing when we are dealing with complicated 
three-dimensional expressions for electric and magnetic fields, and we shall use this principle later on 
to define non-physical potential functions that can describe (via derivatives) complete three-dimensional 
vector functions. 

In the same way, Eq. ( [T5] ) can (and will) be used to describe divergence-free fields with simple 
'vector potentials'. 



1.3 Useful theorems by Gauss and Stokes 

The theorems of Gauss and Stokes are some of the most commonly used transformations in this chapter, 
and therefore we shall take a moment to explain the concepts of them. 



Gauss 's theorem 

Gauss's theorem not only saves us a lot of mathematics but also has a very useful physical interpretation 
when applied to Maxwell's equations. Mathematically speaking, we transform a volume integral over 
the divergence of a vector into a surface integral over the vector itself: 

J dV = ^a-dS. Gauss's theorem (16) 

V 'sources' s 



The surface on the right-hand side of the theorem is the one that suiTounds the volume on the left-hand 
side. If we remember that the divergence of a vector field is equal to its sources, Gauss's theorem tells 
us that: 



- The vector flux through a closed surface equals the sources of flux within the enclosed volume. 

- If there are no sources, the amounts of flux entering and leaving the volume must be equal. 

These statements can be applied directly to Maxwell's equations. Using Eq. (|3]) and applying 
Gauss's theorem, we obtain 

V-EdV = ^E-dS = ^ (17) 

V s 

(Fig. 1), which means that one can calculate the amount of charge in a volume simply by integrating the 
electric flux lines over any closed surface that surrounds the charge, or vice versa. 



Fig. 1: Example of electric flux lines emanating from electric charge in the centre of a sphere 

The same trick can be applied to the source-free magnetic field. Here, we use Eq. Q and obtain 

yV-BdV = ^B-dS = 0. (18) 
V s 



Equation ( [T8| ) gives us the proof of what was already stated earlier: magnetic field lines have no sources 
(V • B = 0), and therefore the magnetic flux lines are always closed and have neither sources nor sinks. 
If magnetic flux lines enter a volume, then they also have to leave that volume (Fig. 2). 



Fig. 2: Example of magnetic flux lines penetrating a sphere 



Stokes's theorem 

Whereas Gauss's theorem is useful for equations involving the divergence of a vector, Stokes's theorem 
offers a similar simplification for equations that contain the curl of a vector. With Stokes's theorem, we 
can transform surface integrals over the curl of a vector into closed line integrals over the vector itself: 

J (Vxa)-dA = ^a-dl. Stokes's theorem (19) 

A C 

One can interpret Stokes's theorem with the help of Fig.[3]as follows: 



- the area integral over the curl of a vector field can be calculated from a line integral along its closed 
borders, or 

- the field lines of a vector field with non-zero curl must be closed contours. 



T JLi T 1 
U J U J 

Fig. 3: Illustration of Stokes's theorem 

The meaning of these statements becomes immediately clear when we apply Stokes's theorem to 
Maxwell's equation ([T]): 

(20) 



J (VxH)-dA = ^H-dl = J 



In the electrostatic case, the time derivative disappears and the area integral over the current density may, 
for instance, be the current flowing in an electric wire as shown in Fig. [4] This means that with a one-line 
manipulation of Maxwell's equations, we have derived Ampere's law, which tells us that every current 
induces a circular magnetic field around itself, whose strength can be be calculated from a simple closed 
line integral along a circular path with the current at its centre. 




->• I ^ H • dl = / Ampere's law (21) 

c 

Fig. 4: Illustration of Ampere's law 

With similar ease, we can derive Faraday's induction law, which is the basis of every electric motor 
and generator. All we have to do is apply Stokes's theorem to Maxwell's equation ([2]): 

(V X E) • dA = ^ E • dl = - ^ y B • dA , (22) 




dt 

and again, after one line, we obtain one of the fundamental laws of electrical engineering. 

Faraday's law tells us that an electric voltage is induced in a loop if the magnetic flux i/A penetrating 
the loop changes over time, as shown in Fig. |5] Alternatively, one can change the flux by moving the 
loop in or out of a static magnetic field. 

I hope that these examples have convinced you that Maxwell's equations are indeed very power- 
ful, and that with a bit of vector analysis we really can derive everything we need for RF engineering 
(although maybe not always in one line . . . ). 

1.4 Displacement current 

Although most people have an idea of what electric and magnetic fields are, the displacement current 
dD/dt is often not so well understood. Since it is vital for the propagation of electromagnetic waves. 
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Faraday's induction law 



(23) 



Fig. 5: Illustration of Faraday's induction 
law 



we shall spend a few lines studying this quantity. We start by deriving and interpreting the continuity 
equation, and then look at a simple practical example. 

We apply the divergence to Maxwell's equation ([T]): 

V-(VxH) = V-J + V— . (24) 
V ' df 

= 




Using Maxwell's equation ([3]l, we have made an association between the 'sources of the displacement 
current' V- (dD/df) and the 'rate of change of electric charge' (d/df)pv. Using the identity ( [T5] ), we 
obtain the continuity equation 

V • J = — , continuity equation (25) 



to which we apply a volume integral and Gauss's theorem ( |T6| ): 

y" V • J dV = ^ J • dS = £ 4 = - ^ y" Pv dV . continuity equation 

V S V 

In this form, the interpretation is very straightforward, and we can state that: 



(26) 



- if the amount of electric charge in a volume is changing over time, a current needs to flow; or, 
more poignantly, electric charges cannot be destroyed. 

Now, it is good to know that electric charges cannot be destroyed, but that does not yet help us 
to understand the displacement current. For this purpose, we go back to Eq. ( |24l ) and this time we do 
not replace the expression for the displacement current. Instead, we apply a volume integral and Gauss's 
theorem and obtain 

fj-dS = Y,In = -^JiydS = -^Jp,.dV, (27) 

S S V 

which we can apply to the simple geometry of a capacitor shown in Fig. [6j which is charged by a static 
current /. 

If we assume a small volume with a surface 5 around one of the capacitor plates, then we can 
directly interpret Eq. ( pTj ): the current /, which enters the volume V on the left, equals the flux integral 
of the displacement current — (d/d?)D, which leaves the volume V on the right. This means that the 
displacement current can be understood as a 'current without charge transport', which in this case can 
only exist because of the rate of change of the electric charge (—{d/dt) J PvdV) on the left capacitor 

V 

plate. 




Fig. 6: Example of a displacement current: charging of a capacitor 



1.5 Boundary conditions 

Before we try to calculate electromagnetic fields in accelerating cavities, we need to understand how 
these fields behave close to material boundaries, for example the electrically conducting walls of a cavity. 
Using Stokes's and Gauss's theorems, we can quickly derive these boundary conditions. 



Field components parallel to a material boundary 

We start with the field components (£"11, //||) parallel to a surface between two materials, as depicted 
in Fig. [7] We define a small surface AA, which is perpendicular to the boundary and encloses a small 
cross-section of the boundary area. Then we integrate Maxwell's equations ([T]l and (|2]) over this area and 
apply Stokes's theorem: 



y^VxH-dA = ^H-dl = y^J-dA + ^yD-dA, (28) 



c 




VxE-dE=<fE-dl=- - / B-dA. (29) 
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Material 2 



Fig. 7: Boundary conditions parallel to a material boundary 



Using Stokes's theorem, the area integrals over A are transformed into line integrals around the 
contour C of the area. If the width d of the area (see Fig. [7]) is now reduced to zero, the calculation of the 
contour integral simplifies to a multiplication of the field components E\\ and H\\ by the path elements 
AZ. The area integrals over D and B vanish and the area integral over the current density J is replaced 
by a surface current, which may flow in the boundary plane between the two materials, times the path 



element AZ. This results in the following boundary conditions: 
//ill —Hn-i = i' , 



conditions for magnetic and 



1^ | 2 ' electric fields parallel to a (30) 

^1 1 = ^1 2 • material boundary 

In the case of a waveguide or an accelerator cavity, we generally assume one of the materials (e.g., 
material 2) to be an ideal electrical conductor, and in that case the electric and magnetic field components 
in this material vanish, so that we obtain 

rr ./ conditions for magnetic and 

Hui = I , 

' electric fields parallel to ideal (31) 

1 1 ~ " ■ electric surfaces 



Field components perpendicular to a material boundary 

In a very similar way, we can derive the boundary conditions for fields (D±, B±) that are perpendicular 
to a boundary surface between two materials. This time, however, we do not define an area but a small 
cylinder with a volume AV around the boundary, as shown in Fig. [8] We form a volume integral from 
Maxwell's equations Q and (|4]) over the volume of the cylinder and apply Gauss's theorem to transform 
the volume integrals into surface integrals: 

Jv-DdV = <j^D-dS = J qydV , (32) 

V S V 

V-BdB = ^B-dS = 0. (33) 
s 



D±i, B±i 




AS2 

Fig. 8: Boundary conditions perpendicular to a material boundary 

In the following step, we reduce the height of the cylinder to zero, so that we end up with two 
surfaces, one on each side of the boundary. And now it becomes clear why we have to start with a 
volume integral. Since the surface element dS is perpendicular to the surface of the cylinder, the 'dot 
product' in the integrals basically reduces the vector fields D and B to the components perpendicular to 
the surface of the cylinder. This means that above equations can now be written as 

conditions for dielectric 
D±i — D±2 = qs , displacement and magnetic 

D , , — R , ^ induction perpendicular to a 

material boundary 

where is a surface charge (measured in units of C/m^) that may exist on the boundary surface 



(34) 
. In die 



case where material 2 is an ideal conductor, we obtain 

conditions for dielectric 
^±1 = , displacement and magnetic ^^5^ 

= . induction perpendicular to 

ideal electric surface 

We note that when the fields are parallel to a boundary surface, the electric and magnetic fields are 
used in the boundary conditions, whereas when they are perpendicular to the boundary surface, we have 
a condition for the dielectric displacement and the magnetic induction. This means that, for instance, the 
tangential electric field £"11 may be smooth across a boundary but there will be a jump in the dielectric 
displacement D|| if there are different relative dielectric constants in the two materials. Similarly, 
the component of the magnetic induction B± perpendicular to a surface may be smooth, whereas the 
magnetic field H± will jump if the two materials have different relative magnetic field constants fJ.^. 



2 Electromagnetic waves 

In this section, we shall derive the general form of the wave equation and then restrict ourselves to 
phenomena that are harmonic in time. Since RF systems mostly deal with sinusoidal waves, we shall 
be able to explain and understand most of the relevant phenomena with this approach. This includes the 
'skin effect', the propagation of energy, RF losses, and acceleration via travelling waves. 



2.1 The wave equation 

We start with the simplification of looking only at homogeneous, isotropic media, meaning we assume 
that the electromagnetic fields 'see' the same material constants (/x, £, k) in all directions. With this 
assumption. Maxwell's equations can be conveniently expressed in terms of only E and H: 

dE 

VxH = JcE + e— , (36) 
ot 

dU 

VxE = -/i^r— , Maxwell's equations (37) 

ot 

V-E = j, (38) 
V-H = 0. (39) 



The curl of Eq. p7) together with Eq. ( [36| ), and the curl of Eq. p6[ together with Eqs. pT\ and ( [38) 
result in the general wave equations for a homogeneous medium 

V^E - V ( V . E) = M .^E + MS , ^^^^ j„ ^ 

H (f homogeneous medium 

dt dt^ 

In the case of waveguides and cavities, we can simplify these equations even further by considering only 
the fields inside the waveguide or cavity, which exist in a non-conducting medium ((C = 0) and a charge- 
free volume (V • £ = 0): 

V^E = U£— E wave equations in a 

non-conducting, charge-free (41) 

y 2 jj _ „ £ jj homogeneous medium 

df2 



2.2 Complex notation for time-harmonic fields 

The already compact wave equations in Eq. ( |4T] ) can be simplified even further by taking into account 
the fact that in RF engineering one usually deals with time-harmonic signals, which are sometimes 
modulated in phase or amplitude. We can therefore introduce the complex notation for electric and 
magnetic fields. We start by assuming a time-harmonic electric field with amplitude Eq and phase (p, 

E{t) =Eo cos {cot + q)) , (42) 

which we can interpret as the real part of a complex expression, 

E{t)=3i{Eoe"Pe"^'] = 3i{Eo cos (cot + q)) +iEo sin {cot + q))} . (43) 

In this form, we can easily separate the harmonic time dependence cot from the phase information (p. 
The phase information can be merged into the amplitude by defining a 'complex amplitude' or 'phasor' 

E = Eoe"P . (44) 

We keep in mind that the real physical fields are obtained as the real part of the complex amplitude times 

£ocos {cot + (p) = {Ee''^'} . (45) 

To simplify our writing, we skip the part with the harmonic time dependence and omit the tilde, 
which means that from now on all field quantities are written as complex amplitudes. In order to convince 
you that this really is a simplification, let us consider what happens to time derivatives when complex 
notation is used: 

-Ee"^' = icoEe"^' . (46) 
dt 

This means that all time derivatives in Maxwell's equations and also in the wave equations can 
simply be replaced by a multiplication by ico, and we are able to do this because the time dependence is 
always harmonic. Only when we have to deal with transient events, such as the switching on of an RF 
amplifier or the sudden arrival of a beam in a cavity, do we have to go back the non-harmonic general 
equations. 

As our first application of the complex notation, we rewrite Maxwell's equations as follows: 

VxH = /weE, (47) 
V X E = -iconH , Maxwell's equations in (48) 

V E = ^ , complex notation (49) 

V-H = 0, (50) 

where the complex dielectric constant e is defined as 

e = e' — ie" = e(l — i—^. complex dielectric constant (51) 

V coe/ 

We note that e is complex only in a conducting medium. We can now proceed to write the general wave 
equations in complex form: 



V^E - V (V • E) = -k^E , 



general complex 
wave equations 



(52) 
(53) 



Here also, we note that the complex wavenumber ^ becomes real in the case of a non-conducting medium: 



1^ = (tp-^e = (O^ne ^1 - ■ complex wavenumber (54) 

Finally, we simplify the wave equations again for the case of a non-conducting, charge-free medium and 
obtain 

V^E = —k^E , complex wave equations in a non- 

(55) 

V^H = —k^H , conducting, charge-free medium 

(56) 

with 

k = CO LIE = . free-space wavenumber (57) 

On the way, we have also introduced a simple definition for the speed of light, c = 1 / y^/l£, in 

Eq. 

2.3 Plane waves 

As an introduction to the theory of electromagnetic waves, we look at a very simple case, that of so-called 
plane waves. We assume again that we are in a homogeneous, isotropic, linear medium and that there are 
no charges or currents, which means that Eqs. (|52]) and ( [53| ) apply. Furthermore — for a plane wave — we 
assume that the field components depend only on one coordinate (e.g., z). The solution of the harmonic 



wave equations (52 1 and (53 1 can then be written as a superposition of two waves 



E,{z)=C,e-T' + C2en^, 

1 (58) 
Hy{z) = -{C,e-T^ + C2en^), 

one of which propagates in the positive and one in the negative z direction. The complex propagation 
constant 7 has a real component a, which describes the damping in a lossy material, and a complex 
component /jS, which describes the propagation of the wave. The relation between the propagation 
constant 7 and the wavenumber k is 

Y=a + i[5 =ik = m^JJie. propagation constant (59) 

We already know that time-harmonic electric and magnetic fields are linked via Maxwell's equa- 
tions, which means that their amplitudes have a certain fixed ratio to each other. This ratio has been 
introduced in Eq. ( [58] ) as the wave impedance Z, the ratio between the electric and magnetic field ampli- 
tudes 

Z = ^ = complex wave impedance (60) 

which becomes real in the absence of lossy material. The wave impedance of free space is given by 

— ss377n. free-space wave impedance (61) 

£0 



2.4 Skin deptli 

When electromagnetic waves encounter a conducting (lossy) material, we have to evaluate the bound- 



ary conditions (see Section |1.5| ), and we find that the wave amplitudes are attenuated suddenly by an 

' case we can assume that 

> 1 , (62) 



attenuation constant a. In the RF case we can assume that 

K 

(oe 



which means that the complex wavenumber ( 54 1 and, obviously, also the complex dielectric constant 



( 5 1 1 are dominated by their imaginary parts, so that we can write 



-ie" = -i— or l^ = —i(0}iK, (63) 
CO 



which is actually equivalent to neglecting the displacement current. Using Eq. ( [59| ), we can then write 
the propagation constant as 



y=a + i^ =ik = m'J = {\+i)J-^ , (64) 

which defines the attenuation constant a. The 'skin depth' is then defined as the distance after which the 
wave amplitudes have been attenuated by a factor 1/e f« 36.8%: 



1 



a y (OjJLK 



skin depth (65) 



Knowing the value of the skin depth is crucial for the design of RF equipment. Let us assume that 
we want to build an accelerating cavity that resonates at 500 MHz. Since high-quality copper is quite 
expensive, we consider the possibility of constructing the cavity out of steel and then copper-plating the 
interior in order to obtain a good quality factor and reduce the losses in the surface. From Eq. ( |65l ), we 
calculate that the skin depth in copper is approximately 3 |am. Depending on how well the copper plating 
is done by the plating company, we can now define the thickness of the copper layer that is needed on 
the inside of the cavity. Typically, around 10-20 times the skin depth is chosen as the plating thickness. 
Figure [9] shows the dependence of the skin depth on the RF frequency. 

Furthermore, the skin depth allows us to calculate the losses in the surface easily. For a wave 
travelling parallel to a conducting surface, one can define a surface resistance by assuming a constant 



current density in a layer of the surface material equivalent to the skin depth, as shown in Fig. 10 



^surf = [^1 • surface resistance (66) 

K-ds 

This value has to be multiplied by l/w to obtain the full RF resistance, where I is the length of the 
conducting wall and w is its width. 

2.5 Energy and transport of energy 

We start this section by presenting Poynting's law, and then explain its components. Poynting's law states 
nothing more than the conservation of electromagnetic energy: 

-^^JwdV = Js-dA + jE-JdV. Poynting's law (67) 
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Fig. 9: Skin depth versus RF frequency 
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Fig. 10: Skin depth and surface resistance 



This equation, read from left to right, states that 'the rate of change of stored energy in a volume 
equals the energy flow out of the volume (through a surface A surrounding the volume) plus the losses 
within the volume (the work performed on charges per unit time)'. In the following lines, we shall see 
that the components of Poynting's law do indeed correspond to what is stated in the previous sentence. 

What is E J? 

In order to understand the expression E • J, we follow |[T| and start with the force acting on a charge 
moving in an electromagnetic field, 

F = g (E + V X B) . Lorentz force (68) 

Multiplying this equation by v and knowing that a • (a x b) = 0, we obtain an expression for the work 
done on a charge per unit time, 

vF = ^vE. (69) 
Assuming N particles per unit volume, we can write 



A^vF = A^^7vE = J-E. 



(70) 



Therefore the expression J • E must be equal to the work done on charges per unit time and unit volume, 
or, in other words, the loss of electromagnetic energy per unit volume. 



The Poynting vector S and the energy density w 

These quantities can be understood by manipulating Maxwell's equations (compare, e.g., We mul- 
tiply Eq. Q by E: 

E-J = E-(VxH)-E-^. (71) 

ot 



Using Eq. ( |D.1[ ), this can be rewritten as 



E-J = H-(VxE)-V-(ExH)-E- — . (72) 

ot 

Using the second of Maxwell's equations ^ and assuming time-invariant jj. and e, we can write 

E-J = -V-(ExH)-^ Qe-D + ^H-B^ . (73) 

Applying a volume integral together with Gauss's theorem ([16]) and rearranging the elements of 
the equation, we end up with 

-I f C-E.D+'-n.B) dV 

"f-lv^^ 2 J Poynting's law (74) 

ExH)-dA+ / E-JdV, 

A Jv 



which can be compared directly with Eq. ( [67] ). On the left-hand side we have the definition of the energy 
density, 

1 1 electric and magnetic energy 

W = Wel + W„,ag = -E • D + -B • H , ^^^.^^ (75) 

and from the right-hand side we obtain the definition of the energy flux density, or the Poynting vector S, 
S = E X H . Poynting vector (76) 

The Poynting vector gives us the direction in which an electromagnetic wave transports energy, 
and from the cross product we understand that this direction is always perpendicular- to the electric and 
magnetic field components. This is consistent with Section [23| where we found that the field components 
(Ex, Hy) of a plane wave (see Eq. ( [SS] )) are perpendicular to the direction of propagation (z). 

In the above derivation, we have used Maxwell's equations in their general form, meaning with 
time derivatives. In the case of the complex notation, the definitions of the energy density and Poynting 
vector have to be modified as follows (for a proof, see ||2| or [^Sjl): 

1 1 electric and magnetic energy ^^^^ 

W = Wel + Wmag= -E-D* + -B-H*, A , ■ T 

*5 4 4 ' density in complex notation 



S = ^ (ExH*) 

2 ^ ' 



complex Poynting vector 



(78) 



3 Electromagnetic waves in waveguides 

In this section, we derive the field components of electromagnetic waves that propagate in waveguides. 
The same principle can then be used to calculate the standing-wave pattern in an accelerating cavity, 
which is nothing more than a superposition of two waves travelling in opposite directions. 

3.1 Classification of modes in waveguides and cavities 

Before we start to solve the wave equation, we need to introduce a classification of the field patterns that 
can be found in waveguides and cavities. 

TMmnp modes, or E,„np modes 

These modes have no magnetic field in the direction of propagation (z) and are therefore often called 
transverse magnetic, or TM, modes. On the other hand, they have an electric field component that is 
parallel to z, hence the equivalent name E modes. 

The indices m, n, p indicate the number of zeros or variations in the three directions of a coordinate 
system. In the case of a waveguide, only the first two indices are used, whereas in the case of a cavity, 
owing to the standing-wave pattern along z, all three are needed for a complete description. In the case 
of a circular waveguide or cavity, the indices indicate the following: 

m, number of full-period variations of the field components in the azimuthal direction. For circularly 

symmetric geometries, E, B oc cos(m<p), sin(m<p). 
n, number of zeros (Xmn) of the axial field component in the radial direction. For circularly symmetric 

geometries, Ez,B^oc J„{xnnr/Rc)- 
p, number of half-period variations of the field components in the longitudinal direction, with E, 

B oc cos{pnz/l), sin{pnz/l)- 

The functions /,„ introduced above are Bessel functions of the first kind and of mth order, and can be 
found in mathematical textbooks. The first three orders are shown in Fig. 11 
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Fig. 11: Bessel functions of the first kind up to order 2 




TEfnnp modes, or H,n„p modes 

Here, there is no electric field in the direction of propagation z, hence the name transverse electric, or 
TE, modes. In analogy to the E modes, H modes have a magnetic field component parallel to z. The 
indices have the same meaning as above. 

TEM modes 

This class of modes has neither an electric nor a magnetic field component in the direction of propagation. 
They can exist between two isolated conductors, for example in a coaxial line. The advantage of TEM 
modes is that waves of any frequency can propagate, whereas TE and TM modes always have a cut-off 
frequency, below which they are damped exponentially (more on this later). However, the disadvantage 
of coaxial lines is that the losses in the two conductors are generally higher than in rectangular or circular 
waveguides. 

3.2 Solution of the wave equation in a cylindrical waveguide 

Instead of trying to find solutions for all six vector components of the electric and magnetic fields, one 
can simplify the problem by using a vector potential A (without any physical meaning) that has only one 
component. One can then quickly derive all six field components from this vector potential. 

It can be shown that only two types of modes can exist in waveguides: TM and TE modes, as 
introduced above. For each mode type, we introduce a vector potential A as follows. Since H and E are 
divergence-free, and since V • (V x a) = 0, we can write 

H™ = VxA™ with E™ = — ^Vx(VxA™), vector potential for TM waves (79) 

(oe 

E™ = V X A™ with H"^^ = —V x (V x A"^^) . vector potential for TE waves (80) 

(OfX 

In both cases the vector potential obeys the wave equation 

y^X = -k^X with k^ = (o^ne, (81) 

which can then be solved for various coordinate systems and has only one vector component, in the 
direction of propagation: 

A=A,e,. (82) 




Fig. 12: Geometry of a circular waveguide 



Circular waveguides 



In a circular waveguide, as shown in Fig. 12 the vector potentials for the TE and TM modes are identical: 



^TM/TE _ cos{m(p)e^''^^^ , vector potential for circular waveguide (83) 



with 



k^ = \Jk'^ —k^. wavenumber in z direction 

Using Eq. (|79]l, we can derive the field components for the TM modes: 



(84) 



Er 



9 



i dH(p 
COS dz 
i dHr 

coe dz 
ik^ 

—A, 
COS 

IdA, 
r d(p 



k k 

■ -C—J'^{kcr)cos{m(p) 

(0£ 

inky 

■ C Jmikcn sm(m(p) 

(oer 

ik^ 

■ C — -Jm(kcr) cos(m(p) 

m 

■ —C—Jm{kcr) sin(m<p) 



-CkcJ'„{kcr) cos(m<p) 



> e 



±ik.z 



field components for TM 
modes in a circular waveguide 



(85) 



Now we can use the boundary conditions to specify the cut-off wavenumber k^. From Section [T3j 
we know that the electric field components parallel to the waveguide surface have to vanish at the surface, 
which means 

a 



Ez{r = a) 







(86) 



The nth zeros of the Bessel functions of order m are tabulated in mathematical textbooks (e.g., Q). 
Using 



271 



c 



(87) 



we can define the cut-off frequency of the waveguide, 

]mn 

(On = C . 



cut-off frequency for TM 
modes in a circular waveguide 



(88) 



The mode that is most commonly used in circular waveguides is the TMoi mode, which has only three 
field components. By inserting m = and n = I into Eq. ( [85] ) and using /o('') = ('')' obtain 
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-C^h{k,r) 
(0£ 



> e 



±ik.z 



■ CkcJi (kcr) 



field components of TMqi 
mode in a circular waveguide 



(89) 



with a cut-off frequency ft)c ~ c x (2.405/a). 

The field pattern of the TMqi mode is shown in Fig. 13 for a mode frequency 15% above the cut- 
off frequency. The distance between the minima or between the maxima of the field corresponds to 0.5 
times the propagation wavelength A,. With decreasing mode frequency. A, = 2n/k, becomes longer, and 
finally becomes infinite when the mode frequency equals the cut-off frequency (Oc- This effect is shown 



in Fig. 14 where the TMqi mode propagates at a frequency just 0.5% above the cut-off frequency. 
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Fig. 13: Field lines of a TMqi mode in a circular waveguide with (O = 1.15(0^. Solid lines, electric 
field lines; dashed lines, magnetic field lines. The brightness of the background is proportional to 
the norm of the field vector: light areas indicate high-field regions of the magnetic field in the left 
plot and of the electric field in the right plot. 
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Fig. 14: Field lines of a TMqi mode in a circular waveguide with (O = l.OOSOc. Solid lines, electric 
field lines; dashed lines, magnetic field lines. The brightness of the background is proportional to 
the norm of the field vector: light areas indicate high-field regions of the magnetic field in the left 
plot and of the electric field in the right plot. 



Rectangular waveguides 

The derivation of the fields in a rectangular waveguide follows the same principle as that used in the 



previous section for circular waveguides. In a rectangular waveguide, as shown in Fig. 15 two different 
vector potentials are needed to describe the TE and TM modes: 




Fig. 15: Geometry of a rectangular waveguide with transverse dimensions a and b 



, ., vector potential for TM waves 

A, = Csin(/c,x) smikyy)e ' =^ j„ ^ rectangular waveguide 



aJ^ = Ccos{kxx)cos{kyy)e 



±ik,z 



where 



k, = ^lB-kl, with kl = kl + k] . 



vector potential for TE waves 
in a rectangular waveguide 



wavenumber in z direction 



(91) 



(92) 



We note that the position of the origin of the coordinate system is linked to the sine and cosine terms in 
Eqs. (90 1 and ( |9T] ). The fields derived from the vector potentials have to fulfil the boundary conditions 
on the waveguide walls. So if, for instance, we were to choose the origin in the centre of the waveguide, 
then the sine and cosine expressions would have to be exchanged to account for the changed symmetries 



with respect to the coordinate axes. Using Eq. ( [79| ) again, we derive the field components for the TM 
modes: 
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: Cky sin(^;cx) cos{kyy) 
- —CkxCOs{kxx) sm{kyy) 
Using the boundary conditions, we can specify the wavenumbers kx and ky-. 

kr = '— and m = 0,1,2,..., 



Ey{x = a) = 
E^lx = a) = 



mn 
a 



Ex{y = b)=0 
E,{y = b)=0 



mi 



and the cut-off frequency for a rectangular waveguide is 

: ckc = cJkl + k^ = cny 



(Or 



m\2 /w\2 

- + 
a J 



and « = 0, 1,2, . . . , 



cut-off frequency for TM 
modes in a rectangular 
waveguide 



(94) 
(95) 

(96) 



The usual convention is to have a> b, and in this case the TEio mode is the mode with the lowest 
cut-off frequency. It is also the only mode that propagates in a relatively large frequency band, from 
/J^''" to 2/7^''°, which is why it is the mode most commonly used in rectangular waveguides. The 
fields of the TE modes can be derived from the TE vector potential using the same procedure. 



3.3 Wave propagation and dispersion relation 



In Figs. 13 and 14 we have seen that the propagation wavelength of a waveguide mode is determined 
by its frequency and by how far the mode frequency is above the cut-off frequency of the waveguide. If 
the propagation wavelength depends on the mode frequency, we can assume that the phase velocity of a 
particular mode also depends on the mode frequency. This relationship is called the dispersion relation 
and, using the definition of the wavenumber in Eq. we can write 



k=k —k= r — = , dispersion relation (97) 

from which we can immediately see that: 

- can be real only if the mode frequency CO is above the cut-off frequency co^ ; 

- for w < ft)c, the mode cannot propagate and the fields are exponentially damped. 

We also have a definition of the phase velocity, which is the speed at which the maxima and minima of 
the field patterns move along the waveguide: 



phase velocity 



(98) 



This is not to be confused with the speed with which the wave actually propagates in the waveguide. The 



dispersion relation is usually plotted in the form of a 'Brillouin diagram', as shown in Fig. 16 




Fig. 16: Dispersion relation in a waveguide. The dotted line shows the case Vph = c. 



The slope of the dispersion relation is the group velocity 



Vgr = — , group velocity (99) 

&kz 

which gives the velocity with which a signal or energy is transported along the waveguide. From Fig.[T6j 
we can conclude that: 

- Each frequency has a certain phase velocity and group velocity, which means that signals with a 
frequency bandwidth will become deformed while travelling along a waveguide. With the help of 
the dispersion relation, we can easily quantify how much deformation will occur. 

- The phase velocity Vph is always larger than the velocity of light c, and at cut-off (CO = (0^) it even 
becomes infinite {k^ = and Vph — )■ °°). 

- For acceleration, one needs synchronism between the phase velocity (the speed of the field pattern) 
and the velocity of the particles, which implies that acceleration in waveguides is impossible. 

- Information and therefore energy travel at the group velocity, which is always slower than the 
speed of light. 



3.4 Attenuation of waves (power loss method) 

Up to this point, we have assumed perfect electrical conductors as the boundaries of our waveguides. Real 
waveguides and cavities have a certain resistance, and the fields therefore penetrate into the conductors, 
which significantly complicates the solution of the wave equation. However, we have seen in Section 



2.4 that the skin depth in metals is very much smaller than the RF wavelength. This means that we 



can reasonably assume that the field patterns in a waveguide with ideal boundaries and in a waveguide 
with resistive metal boundaries will be practically identical (of course, only for good conductors such as 
copper or aluminium). In order to calculate the attenuation of waves, we can therefore use the fields of 
a waveguide with ideal electrical boundaries. From the magnetic field, we calculate the induced current 
in the waveguide walls, and then apply the resistance of the real material to calculate the losses and then 
the damping of the wave. This principle is called the power loss method and is a simplified method for 
calculating RF losses on the surfaces of good conductors. 

We start by defining the power that is lost per unit length along the longitudinal axis of the wave- 
guide. 



P' = ~. (100) 
dz 

From 

E,H oce^"' Poce-2az^ (-jOl) 

we immediately obtain 

dP 

P' = -— = 2aP (102) 
dz 

and thus the definition of the attenuation constant 



P' 

a = — . attenuation constant (103) 

2P 

In the next steps, we need to derive expressions for the power P transported through the waveguide, 
and the power loss per unit length P'. Using the field components of the TMqi mode given in Eq. ([89]) 
and the definition of the complex Poynting vector in Eq. (78 1, we obtain 



a 27C ^ ^ T T 

,2 72/ 



f = l/(ExHn.dA= i//E,H;.d.dy^ , (104) 







where we have used 



/ 



2 



J^{kcr)rdr='^J^^{k,a). (105) 


In order to calculate the losses on the waveguide surface, we first need to know the surface currents 



that flow within the skin depth. For this purpose, we make use of Ampere's law, as shown in Fig. 17 



^H-Jl = / = ^ J-(5sdl). 

c c 



Ampere's law 



(106) 




r-Htp = ; aAq) 



k=Kai 



Fig. 17: Ampere's law applied to calculate the surface currents in a circular waveguide 



Since the magnetic field has only an azimuthal component, we obtain 

H^{r = a,z)= CkJi {k,a)e-"'^' = J,{z)5s . 



(107) 



The power density (in W/m ) in the waveguide wall is given by 

1 1 d^P 

Pv = -E • r = —7,/; = , power density (108) 

2 2k {dr){rd(p)[dz) 



from which we can write an expression for the power loss per unit length. Together with Eq. ( |107[ ), we 
obtain 

I dp f f , , naC\lj}iha) . ... 

P = ^ = I I p^rdrdm = ^^^-^ — - , power loss per unit length (109) 

OZ J J KOs 

a 

where we have used the fact that 5s ^ a to simplify the evaluation of the integral. Now we insert 
Eqs. ( |104| ) and ( |109| ) into Eq. ( |103| ) and obtain an expression for the attenuation of a TMqi mode in a 
circular waveguide, 

P' /?surf attenuation of TMol mode in 

(^ = ^7^ = = . . , (110) 

2P ZQay'\-{fJfY circular waveguide 



In the expression above, we have used the definition of the surface resistance given in Eq. ( [66| ) and the 
definition of the free-space wave impedance Zq given in Eq. ( [61] ). 

As an example, we have plotted the attenuation constant for an aluminium waveguide in Fig. [18} 
where we can see that for this type of waveguide: 

- Large-diameter waveguides result in smaller losses, which means that a cost optimum has to be 
found between the cost of the waveguide, its space requirements, and the losses. 

- The minimum losses occur when the operating frequency of the TMqi mode is a factor of \/3 
above the cut-off frequency (try to prove this!). 
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Fig. 18: Attenuation of the TMoi mode in a circular aluminium waveguide for several different radii: 
bottom to top, 0.5 m, 0.4 m, 0.3 m, 0.2 m 



4 Accelerating cavities 
4.1 Travelling- wave cavities 

In order to accelerate particles in a 'waveguide-like' structure, the phase velocity in the structure needs 



to be slowed down, which can be achieved by putting some 'obstacles' into the waveguide. In Fig. [19 
we see a simple example of a disc-loaded waveguide. 
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Fig. 19: Geometry of a simple travelling-wave structure 



The dispersion relation for such a structure is derived in, for example, (5) as 
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2.405c 



dispersion relation of 
disc-loaded circular 
waveguide 



where 
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4fl3 



3nJ^{2A05)b^L 



1 and a : 



2.405 



(111) 



(112) 



Plotting Eq. ( |111[ ) gives us the Brioullin diagram in Fig. 20 where we can see that we now obtain phase 
velocities that are equal to or even below the speed of light. We can also understand why the 27r/3 mode 
is often used for acceleration in electron accelerators, because for this mode (in this example) the phase 
velocity is just equal to the speed of light. It should be noted that with different geometries, it is possible 
to operate with different modes and also at velocities Vph < c. When a structure operates in the 27r/3 



mode, this means that the RF phase shifts by 27r/3 per cell, or, in other words, one RF period extends 
over three cells. 




Fig. 20: Dispersion diagram for a disc-loaded travelling-wave structure. Here, the chosen operating 
point is Vph = c and k- = 2n/3L. 



By attaching an input and an output coupler to the outermost cells of the structure, we obtain a 
usable accelerating structure. Since the particles gain energy in every cell, the electromagnetic wave 
becomes increasingly damped along the structure. It is then extracted via the output coupler and dumped 
in an RF load. If one is interested in obtaining the maximum possible accelerating gradient in each cell, 
then one can counteract the decreasing fields by changing the bore radius from cell to cell. The idea is 
to slow down the group velocity from cell to cell and obtain a 'constant-gradient' structure, rather than a 
'constant-impedance' structure where the bore radius is kept constant. Other optimizations, for example 
for maximum efficiency, are also possible. 



4.2 Standing-wave cavities 

One obtains a cylindrical standing-wave structure by simply closing both ends of a circular waveguide 
with electric walls. This yields multiple reflections on the end walls until a standing-wave pattern is es- 
tablished. Owing to the additional boundary conditions in the longitudinal direction, we obtain another 
'restriction' on the existence of electromagnetic modes in the structure. Whereas a longitudinally open 
travelling-wave structure allows all frequencies and all cell-to-cell phase variations on the dispersion 
curve, now only certain 'loss-free' modes (still assuming perfectly conducting walls) with discrete fre- 
quencies and discrete phase changes can exist in a cavity. If RF power is fed in at a different frequency, 
then the fields excited are damped exponentially, similarly to the modes below the cut-off frequency of a 
waveguide. 

The corresponding dispersion relation for a standing-wave cavity can again be found in textbooks 
(see |[5j and also Q). However, it is necessary to pay attention to whether the structure under consider- 
ation has magnetic or electric cell-to-cell coupling and what kind of end cell is assumed in the analysis. 
The most common form of the dispersion relation is derived from a coupled-circuit model with N + I 
cells. Usually the model has half-cell terminations on both ends of the chain, representing the behaviour 
of an infinite chain of electrically coupled resonators (compare the original paper by Nagle et al. ||7|): 
^ dispersion relation for 

^" = /. , n = 0,l,...,N. half-cell-terminated (1 13) 

^ standing-wave structure 

Assuming an odd number of cells, coo is the frequency of the n/2 mode and of an uncoupled single cell; 
k is the cell-to-cell coupling constant, and nn/N is the phase shift from cell to cell. For k<^\, which is 
usually fulfilled, the coupling constant is given by 



^iz mode mode 
(Oo 



coupling constant 



(114) 



Two characteristics of the dispersion curve are worth noting: 

The total width of the frequency band of the mode, (Ojc mode — ^^o mode> is independent of the number 
of cells, which means that we can determine the cell-to-cell coupling constant by measuring the 
complete structure (but this is only true if all coupling constants are equal). 
For electric coupling, the mode has the lowest frequency and the 7i mode has the highest. In 
the case of magnetic coupling, this behaviour is reversed, and one can find the corresponding 



dispersion curve by changing the sign before the coupling constant in Eq. ( 113 1 



In Fig. |2T] we plot the dispersion curve for a seven-cell (half-cell-terminated) magnetically coupled 
structure according to Eq. ( |113| ). 




Fig. 21: Dispersion diagram for a standing-wave structure with seven magnetically coupled cells 



In practice, one usually has cavities with full-cell termination, and in this case one has to detune 
the frequences of the end cells to obtain a flat field distribution in the cavity [ 8] . In this case it is possible 
to have a flat field distribution for either the mode or the n mode but not for both at the same time, 
because the end cells have to be detuned by different amounts in the two cases ||9|. 

4.3 Standing wave versus travelling wave 

The principal difference between the two types of cavity is in how and how fast the cavities are filled with 
RF power. Travelling- wave structures are filled 'in space', which means that, basically, cell after cell is 
filled with power. For the following estimations, we assume a frequency in the range of hundreds of 
megahertz. The filling of a travelling-wave structure typically takes place with a speed of approximately 
1-3% of the speed of light and results in total filling times in the submicrosecond range. Standing-wave 
structures, on the other hand, are filled 'in time': the electromagnetic waves are reflected at the end 
walls of the cavity and slowly build up a standing-wave pattern of the desired amplitude. For normal- 
conducting cavities, the time required for this process is typically in the range of tens of microseconds. 
For superconducting cavities, the filling time can easily go into the millisecond range (depending on the 
required field level, the accelerated current, and the cavity parameters). This means that for applications 
that require very short beam pulses (< 1 ^s), travelling- wave structures are much more power-efficient. 
For longer pulses (> « x 10 ps), both types of structures can be optimized to achieve similar efficiencies 
and costs. 



Since one can have extremely short RF pulses in a travelling-wave structure, one can obtain much 
higher peak fields than in a standing-wave structure. This is demonstrated by the accelerating structures 
for CLIC |10|, which have reached values of approximately 100 MV/m (limited by electrical break- 
down), whereas the design gradient for the superconducting (standing-wave) cavities for the ILC 111] is 
just slightly above 30 MV/m (this value is generally limited by field emission and by quenches caused 
by the peak magnetic field). 

Travelling-wave structures can, theoretically, be designed for non-relativistic particles. In existing 
accelerators, however, they are mostly used for relativistic particles. Low-beta acceleration is typically 
performed with standing-wave cavities. 

Because of the lack of an obvious criterion (other than the pulse length or the particle velocity), 
an optimization and costing exercise has to be performed for each specific application in order to decide 
which structure is more efficient. Two excellent papers 1 12] 13 1 in which this exercise is performed can 
be used as references. 



4.4 The pillbox cavity 

In this chapter, we shall analyse only the simplest TM-mode cavity, the so-called pillbox cavity. A 
selection of cavities using other mode types is described in p4[ . 

Resonating cavities can be represented conveniently by a lumped-element circuit consisting of an 
inductor (for storage of magnetic energy) and a capacitor (for storage of electric energy). Looking at 



Fig. 22 one can easily imagine how the lumped circuit can be transformed into a cavity. 





Fig. 22: Transition from a lumped resonating circuit to a resonating cavity 



The pillbox cavity is nothing more than an empty cylinder with a conducting inner surface. The 
simplest mode in this cavity is the TMqio mode, which has zero full-period variations in the azimuthal 
direction (9), one 'zero' of the axial field component in the radial direction (r), and zero half-period 
variations in the longitudinal (z.) direction. We can derive the general field equations by using the vector 
potential for a circular waveguide given in Eq. ([83]) and simply superimposing two waves, one propagat- 
ing in the positive z direction and one in the negative z direction: 



^TM/TE ^ cJn,{krr) COs(m<p) (e"'*^^^ + 
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vector potential for travelling 
waves in the positive and 
negative z directions 



(115) 



Using Eq. (79 1, we derive the TM field components 
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TM modes in a pillbox cavity (116) 



In the case of standing-wave cavities, the term 'cut-off frequency does not really make sense, 
so we have replaced the symbol kc by k^, indicating that we have a radial dependence of the axial field 
component, which can also be interpreted as a radial wavenumber. 




Fig. 23: Pillbox cavity 



In the next step, we apply the boundary conditions for a pillbox with radius a and length L as 



shown in Fig. 23 We obtain 



Er{z = 0/L),E<p{z = 0/L)=0 



E(p{r = a)., E^{r = a), Hr{r = a) =0 =^ k,- 



Jn 



(117) 
(118) 



In the case of the circular waveguide, the transverse boundary condition made a discrete quantity 
out of (which we now call k,- in the above equations), and thus defined the cut-off frequency. Now, 
with the second boundary in the z direction, we obtain a discrete solution for L also. The two boundary 
conditions together result in a discrete set of frequencies (the dispersion relation) for our pillbox cavity: 
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dispersion relation for TM 
modes in a pillbox cavity 



(119) 



We note that the dispersion relation of a single-cell cavity as given above is different from the 
dispersion relation that can be derived for a multicell cavity, as in the case of Eq. ( |1 13 1. The latter is de- 
rived from a model of equivalent lumped circuits, each representing a cell resonating in the TMqio mode 
and coupled to its neighbours in order to model the behaviour of a multicell cavity, whereas Eq. ( |1 19| ) is 
directly derived from Maxwell's equations and describes a field pattern that is created by the boundary 
conditions of our pillbox. 

The TM mode with the lowest frequency is the TMqio mode, with a frequency 



TM _ 2.405c 



and its field components are 
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Figure 24 shows the field pattern of the TMqio niode, simulated by Superfish 
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Fig. 24: Field pattern of the TMoio mode in a pillbox cavity 



4.5 Basic cavity parameters 

In order to characterize and optimize cavities, we need some commonly used figures of merit, which 
we shall define here in general terms and then apply to our simple pillbox cavity. In the following, we 
assume that we are dealing with an axially symmetric cavity resonating in the TMoio mode. 

4.5.1 Energy gain in a cavity 

For particles traversing a cavity on axis, the electric field generally has the following form: 

E^{r = 0,z,t) = E{0,z)cos{(Ot + (p) , 
which we can use to calculate the energy gain of a particle when it traverses the cavity, 

L/2 

AW = q J E{0,z) cos{cot + (jo) energy gain in a cavity 

L/2 (Panofsky equation) 

= qVoTcos(p = qEoTLcos(p , 
where the cavity voltage is given by 



(122) 



(123) 



L/2 

yo= J E(0,z)dz = Eo, cavity voltage (124) 

-L/2 

and the 'difficult mathematics' has been lumped into the so-called transit time factor 

L/2 L/2 

/ £(0,z)cos(cof(z))dz / E{0,z)sm{(ot{z))dz 

-r -L/2 

T = tan . 

, , . . transit time factor (125) 

/ E{0,z)dz J E{0,z)dz 

-L/2 -L/2 

V ' 

=0 if E{0,z) is symmetric about z=0 

This takes into account the fact that the RF electric field changes during the passage of the particles. It 
gives the ratio between the energy gained in an RF field and in a DC field and is therefore always less 
than 1. We note that the Panofsky equation takes account of the changing velocity of the particles when 
they cross the accelerating gap. This makes the integrals in the above equations difficult to evaluate. 
Assuming that the velocity change of the beam particles during their passage is small, however, one can 
say that 

mf=i(o- = ^, (126) 

V pA 

which changes the expression for the transit time factor to (assuming that E{0,z} is symmetric about 
z = 0) 

L/2 

/ £(0,z)cos(2;rz/j3A)dz 
^ -L/2 transit time factor for small 

(Iz/ ) 



L/2 velocity changes 

/ E{0,z)dz 

-L/2 

The accelerating voltage Vacc is the voltage that the particle 'sees' when crossing the cavity and 
should not be confused with the cavity voltage Vq. We thus define 



Vaec =VoT = EqLT . accelerating voltage ( 1 28) 



4.5.2 Shunt impedance 

The shunt impedance tells us how much voltage a cavity will provide when a certain amount of power 
is dissipated in the cavity walls. This is one of the parameters to be maximized in cavity design, since a 
large shunt impedance reduces the power consumption of an RF cavity. The general definition is 

Vq shunt impedance (linac 

^ ~ ' definition) 



The benefit of a high shunt impedance can easily be diminished by having a small transit time factor, be- 
cause in this case the cavity voltage cannot be used efficiently to transfer energy to the beam. Therefore 
one usually tries to optimize both the shunt impedance and the transit time factor, which explains the 



definition of the effective shunt impedance 



R 



(VoT) 



effective shunt impedance 



(130) 



When comparing multicell structures operating at different frequencies, one is interested less in 
the efficiency per cell (because the cell size depends on, for instance, the frequency chosen) than in the 
efficiency per unit length of the accelerating structure. For this reason, we define 



L 



Pd/L 



sliunt impedance per unit 
length 



(131) 



and 



R 
L 



{EpTf 
Pi/L 



effective shunt impedance per 
unit length 



(132) 



4.5.3 'Linac ' and 'circuit' definitions of shunt impedance 

It turns out that different communities of accelerator experts use different definitions of the shunt impedan- 
ce. Linac experts usually use the definitions presented above, whereas the people who deal with circular 
machines generally use a definition that is derived from the lumped-circuit definition of a resonator (see 
Section |477] ). In that definition, all shunt impedances are exactly half as large, following 

Vn shunt impedance (circuit 

2P^ definition) 

So, before you discuss shunt impedances with anyone, make sure that you are using the same definition. 
In order to mark the difference clearly, we use R^ in this text to identify when the circuit definition is 
being used. 

4.5.4 3 dB bandwidth and quality factor 

The quality factor Q describes the bandwidth of a resonator and is defined as the ratio of the reactive 
power (stored energy) to the real power that is lost in the cavity walls: 

^ ^ ^ WW quality factor (134) 

Aft) Pd 

If a resonator were built with ideal electrical walls (zero electrical resistance), the resonance curve 
would be a delta function at the resonance frequency. So, the bandwidth Aco would be zero and the quality 
factor would be infinite. In reality, even superconducting cavities have a certain surface resistance. 



which is why all our cavities have a certain bandwidth and a finite quality factor. Figure 25 shows a 
typical resonance curve measured with a network analyser. In a measurement of this kind, two antennas 
penetrate the cavity. The first antenna sends an RF signal with a frequency sweep, and the second picks 
up the field level in the cavity. As a result, we obtain a plot of the field level versus frequency. The 
bandwidth is defined as the frequency width of the resonance curve, measured as the distance between 
the points where the field level has dropped by 50% (or —3 dB), as shown in Fig. [25] 

Together with the shunt impedance, one can define another figure of merit, {R/Q), which is used 
to maximize the energy gain in a cavity for a given stored energy: 




{R/Q) (135) 



QJ coW 

{R/Q) is independent of the surface losses of the cavity and can therefore be used to quaUfy the geometry 
of an accelerating cavity. 

4.5.5 Filling time of a cavity 



This section is a short extract from 1 14 1, which can be consulted for more details. The dissipated power 



in a cavity must be equal to the rate of change of the stored energy: 

Qq 

The solution of the above equation can be written as 



W{t) = WQe-^*l\ (137) 
which describes an exponential decay of the stored energy with a 'filling time constant' T, where 



Oct 

T = — - . filling time constant (138) 

OX) 

If the cavity is equipped with a power coupler, we have to consider the 'loaded Q' (which will be derived 
later), and the filling time constant changes to 

2Q\ filling time constant for a 

^' (Do loaded cavity 



In the above definition, the electric field decays exponentially with a time constant 1/t, whereas the 



stored energy decays with a time constant 2/t. Be aware that you can often find textbook definitions of 
the filling time constant where the stored energy decays with a time constant 1 /t. 



4.6 Basic cavity parameters for a pillbox cavity 

As a small exercise, in this section we calculate the cavity parameters that were defined in the previous 
section for a pillbox cavity of length L and radius a. Since the TMqio mode has no z dependence, we can 



simplify the expression for the transit time factor ( 127 1 to 



L/2 

/ £(0,z)cos(27rz/j3A)dz 
-L/2 _ sin(7rL/j3A) 

I EiO,z)dz 

-L/2 



transit time factor of a pillbox 
for small velocity changes 



(140) 



In the case of relativistic particles (j3 « 1) and a cavity length L = A/2, which is often chosen because 
the cavity can then be cascaded into a multicell structure, we obtain 



T = =0.64. 
71 



transit time factor of a pillbox 
for relativistic particles 



(141) 



With real cavities, one usually tries to increase the transit time factor by shortening the accelerating gap. 
This can be done by introducing nose cones on the cavity walls, as shown in Fig. |22] 

We use the power loss method again to calculate the quality factor of our pillbox cavity. To evaluate 
Eq. ( |134 i, we need the stored energy and the power lost in the cavity walls. For the stored energy, we 
obtain 

W=Wel+Wn,ag = 2Wel = 2 J^E-T>*dV. 



(142) 



With 



we obtain 



EqJq 



Mr 



a in L/l 

^ = jJJ J Elj'J^-^ydrdcpdz=lEieonLa^{M) 



(143) 



(144) 



-L/2 

To calculate the dissipated power, we integrate Eq. ( |108| ) over a volume that consists of the inner surface 
of the pillbox times the skin depth: 



L/2 



Pd=;^ / JzJ*, 2nadz+- [ JrJ; Hl{r,z = 0)2Krdr 

E^TlRfMYfCl 2/ . 



(1/4)2 



^0 



-Jt{joi){a + L), 



where we have made use of 



701?- 



Putting everything together, we obtain 

(OW 



Zn(0 La 1 La 



2/?surf L + a 5. L + a 



(145) 



(146) 



(147) 



(148) 



As we can see, the quality factor is a function of the material constants K and /i (which ai^e contained 
in Ps), the frequency, and the geometry of the cavity. We also note that for the same cavity shape, the 
quality factor increases with the frequency in proportion to ^/Q). 

The accelerating voltage in a pillbox cavity is given by 



VoT = EoLT = EoL 



sm{KL/pX) 



accelerating voltage in pillbox (149) 



and is obviously a strong function of the transit time factor. It therefore depends on the gap length L and 
the speed of the particles j3. Owing to their high development costs, superconducting cavities are often 
used over large velocity ranges without changing their cell length, and this results in a velocity-dependent 
acceleration efficiency. Figure 26 shows {R/Q) °^ (Vbr)^ as a function of particle velocity for a five-cell 
superconducting cavity whose geometric cell length corresponds to a particle speed of j8 = 0.65. 



C3 
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'0.50 0.55 0.60 0.65 0.70 0.75 0.80 

Fig. 26: Dependence of {R/Q) on particle velocity for a five-cell superconducting cavity with a 
geometric )3 of 0.65 and a frequency of 704.4 MHz. Upper curve, % mode (see also pS)). 



If the cavity is used over too large a velocity range, one may find areas where the passband mode 
that is closest to the n mode (here, the 4/5?! mode) has a higher acceleration efficiency than the acceler- 



ating mode. These areas are highlighted in Fig. 26 and should be avoided when one is designing a linac. 



One should also be aware that the {R/Q) of the HOMs is highly dependent on the particle velocity. 

Using the expressions for the accelerating voltage VqT (Eq. ( |149| )) and the dissipated power 
(Eq. ( |146| )), we also obtain an analytical expression for the effective shunt impedance. 



R 



(VoT) 



sin(7rL/j8A) 



Pd 7rPsurf-/f(joi) TiL/jiX a{a+L)' 



effective shunt impedance of a 
pillbox 



(150) 



Finally, we calculate the frequency and {R/Q) using Eqs. ( |120| ), ( |150| ), and ( |148| ): 



TM _ 2.405c 



pillbox frequency 



(151) 



2c sin(7rL/j8A) L 



piUbox (R/Q) 



(152) 



As stated before, {R/Q) is indeed independent of any material parameters. However, it does depend on 
the geometry of the cavity and the transit time factor. 

4.7 A cavity as a lumped circuit 

In the field of RF technology, it is common practice to describe the behaviour of cavities, RF transmission 
lines, and couplers with equivalent lumped circuits. In this chapter, we shall present only the treatment 
of a cavity and a coupler, so that one can understand how to get power into a cavity. Descriptions of 
the transmission of RF power and the associated theory of RF transmission lines can be found in many 
textbooks on RF and microwave engineering. We start with the description of a cavity by a parallel LCR 



circuit as depicted in Fig. 27 




Fig. 27: Lumped-circuit equivalent of a resonant cavity 

You may remember that the admittance of a parallel circuit is calculated by adding up the admit- 
tances of the individual elements, which means that we can write the cavity impedance as 



1 



i(0C + l/i(0L+\/R'' 



lumped-circuit cavity 
impedance 



(153) 



At resonance (ft) = (Oq), the imaginary parts cancel each other and the cavity impedance becomes real, 
which means that 



X = (OqL ■■ 



(OoC 



that the resonance frequency is given by 



COo 



1 



and that the power lost in the resonator is given by 



lumped circuit at resonance ( 1 54) 



lumped-circuit resonance 
frequency 



(155) 



1 (VoT) 



2 



lumped-circuit dissipated 
power 



(156) 



The stored energy can be written as 



1 1 (V T)'^ 

W = -C{VoTf = - ^^-V- ' lumped-circuit stored energy (157) 

2 2 (OqL 

and from this we obtain an expression for the quality factor, 
W 

Qq = (Oq — = coqCR'^ = . lumped-circuit quality factor (158) 

Pd (OoL 

Our goal is to relate the lumped elements to the cavity characteristics, and for this purpose we 
multiply Eq. ( |157[ ) by ft) and, together with Eq. ( |154[ ), we obtain 



1 [l (VoT)^ fR''\ \ (R 



(OqC VC IcoqW \QJ 2\Q 



(159) 



From this, we can understand the difference between the 'circuit ohm' and the 'linac ohm', and it also 
provides a lumped-circuit description of a cavity, as summarized in Table [T] As we can see, three quan- 
tities are sufficient to describe a resonator. Instead of using R, L, and C, one can also use the parameters 
COq, Qo, and (R/Q) to completely characterize an RE cavity, as in Table |2] 

Table 1: Lumped-circuit elements of a cavity 



Lumped circuit Field description 



2 



(Oo{R/Q) 
1 /R 



2(0o \Q 



Table 2: Three characteristic quantities of a cavity 

Lumped circuit Field description 

1 2.405c 
(Oq = —j= (pillbox) 



ILC a 

R"" (OqW 

Qo = cooCR' = — Qo = ^ 

(OqL Pd 

R'\ [I 1 /R\ fR\ {VoTf 



4.8 Getting power into a cavity: couplers 

In this section, we shall extend the circuit model to include the power coupler and also extend our basic 
equations to describe the process of coupling power into a cavity. There are two basic types of couplers 



that are used in standing- wave cavities: 



AntennaAoop couplers: here, the coupler is usually some kind of coaxial line, with the outer 
conductor connected to the cavity wall and the inner conductor either penetrating into the cavity 
volume or connected in a loop to the inner surface of the cavity (Fig. 28 1. 

Iris couplers: here, the fields in a waveguide are coupled to the cavity fields via an opening that 
connects the waveguide to the cavity. 



Fig. 28: Example of an antenna-type coupler (left) and a loop-type coupler (right) 



When designing a coupler, one has to keep in mind the principle of reciprocity: the coupler has to 
produce a field pattern in the area of the coupling port that is very similar to the field pattern of the mode 



that will be excited in the cavity. Looking at Fig. 28 one can imagine that an antenna-type coupler would 
be very effective on the end walls of our pillbox, where it would couple electrically to the axial electric 
field lines. On the cylindrical surface of the pillbox, a loop coupler would be a better choice, with the 
loop oriented such that the azimuthal magnetic field penetrates the loop. 



Figure 29 shows an example of a 'tuner-adjustable (waveguide) coupler' (TaCo) 1 16 1, as used for 



the Linac4 1 17 1 cavities at CERN. In this case a short-circuited rectangular waveguide is coupled to a 
standing-wave cavity via a racetrack-shaped coupling iris. The coupling factor (more on this later) here 
is a function of the position of the short circuit (left side), the height of the racetrack-shaped coupling 
channel between the cavity and the waveguide (on the top), the size of the coupling slot, and the position 
of a stub tuner, which is used to fine-tune the coupling. 
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Fig. 29: Waveguide coupler connected to a Linac4 cavity 



In the ideal case, the power coupler is matched to the (beam-)loaded cavity, which means that 
there is no reflected power returning from the cavity towards the RF power source. Here, 'matched' 
means that the coupler acts like an ideal transformer that transforms the impedance of the cavity into 
the impedance Zq of the attached waveguide. To keep things simple, let us assume that the RF generator 



is also matched to Zq so that we can estabhsh a lumped-element circuit as shown in Fig. 30 
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Matched generator Waveguide Coupler " Cavity 
Fig. 30: Lumped-element circuit for RF power source, waveguide, power coupler, and cavity 



Considering the coupler as a transformer, we can write that 

, ^ VoT 



VqT = nVgen 



^gen 

n 



n Z^, 



which means that the cavity impedance 

1 



mC+X/mL+X/R" 



cavity impedance 



is transformed into 



X 



2 /DC ' 



(160) 



(161) 



cavity + coupler impedance (162) 



which is the impedance 'seen' from the waveguide. The stored energy in the resonator, expressed in 
lumped-circuit values, becomes 



C 



C 



and the dissipated power can be written as 



stored energy 



(163) 



Pa = = n 



2 y2 

2 gen 



dissipated power 



(164) 



Now we can define the quality factor of the unloaded cavity with lumped-circuit elements: 

= = coqR'^C . unloaded Q (165) 



When the generator is switched off, not only will the stored energy in the cavity be dissipated in the 
cavity walls, but a power Pex will also leak out through the power coupler, where 



Using /"ex. one can define the quality factor of the external load. The external Q is thus defined as 



2ex = — - = n^O}oZoC . external Q (167) 



4.8.1 Undriven cavity 

In order to understand the power balance and matching for a driven cavity with beam, we start with a 
simple case, assuming that the RF is switched off and that there is no beam in the cavity. The power 
balance is then 

power balance of undriven 

Ptot = Pd+Pex, eavity ^^^^^ 

with which we can define the so-called 'loaded Q of the ensemble of cavity and coupler by 



— = — + — . loaded Q (169) 

Q\ Gex Qo 

The coupling between the cavity and the waveguide is described by the coupling factor j8 , where 

i8 = ^ = 1^ = ^ . coupling factor (170) 

Gex « Zo 

Optimum power transfer between the cavity (+ coupler) and the waveguide takes place when the impedance 
at the coupler input equals the waveguide impedance at the resonance frequency of the cavity. We know 
that the cavity impedance becomes real at resonance, which means that 

Zc = = n^Z^ = n^Zo => j8 = 1 . (171) 

It is important to keep in mind that the 'matching condition' j8 = 1 is only valid for a cavity without 
beam. 

4.8.2 RF on, beam on 

Once we take the beam loading into account, the power needed in the cavity increases and will yield a 
different value for the coupling factor j8 at the point of optimum power transfer. A simple way to intro- 
duce the beam is to treat it as an additional loss in the cavity, which can be added to the power dissipated 
in the cavity walls: 

dissipated power + beam , , ^ . , 
Pdb = Pd + /l. po^er ^1^2^ 

As in the case without beam, maximum power transfer to the cavity is achieved when the input impedance 
of the coupler equals the impedance of the waveguide. This condition yields zero reflection and also 
implies that the power needed in the cavity, /^d (for losses and beam), has to be equal to /"ex as defined 
in Eq. ( |166 1. This means that 



where we have introduced a quaUty factor 2ob for the cavity plus beam. For the matched condition, we 
therefore obtain a couphng factor of 



i8 = l + 



and the following quality factors: 

(OqW Qo 



2ex — Qob 



Pb+Pd 

Qo 

l+j8 2 + 



l+Pb/Pi P 
Qo 



matched coupling factor with 
beam 



external Q with beam 



loaded Q with beam 



(174) 



(175) 



(176) 



In the case of a superconducting cavity, one can generally assume that P\, » Pd, which means that the 
coupling factor for the matched condition can be written as 



Using 



^ Pd Pd 



^ =4eamVorcOS^s 



matched coupling factor for 
SC cavity + beam 



(177) 



(178) 



we can write a simple expression for calculating the loaded and external Q values for a superconducting 
cavity as follows: 



Qi-Q. 



Qo 



VoT 



Pbsam/Pd (^/2)4eamCOS(/)s 

The results in this paragraph are summarized in Table [3] 



Q\/ex for SC cavity 



(179) 



Table 3: Definitions of Q values and coupling factors for driven and undriven cavities 



General 



Matched case 



Undriven cavity 



Driven cavity 
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Gex 


^Qo 
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' ex 


_ Qo 






" Gex 



Pex _ Qo_ 
Pd ~ Gex 

Gex = Qo 



1 



a 

^j8 = l 

(OqW 



Qx 



Qo 



Qo 

^ex _ Qob 
Pdb Gex 



l^j8 = l + 



Gex = Qob 

Qx 



_ (OqW 
~ Pd+Pb 

Qo 



2 + Pb/Pd 



4.9 'Matching' a cavity 

In the last section, it was claimed that part of an electromagnetic wave is reflected when it 'sees' a change 
in impedance during its propagation. In fact, the whole purpose of the power coupler was to transform 
the impedance of the waveguide into the impedance of the cavity. In this last section, we shall see why 



this is so. For this purpose, we look at a transmission line as shown in Fig. 31 
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Zo 

Fig. 31: Voltages and currents along a transmission line 



This transmission line is representative of a waveguide, a coaxial line, or any other kind of trans- 
port geometry used to guide electromagnetic waves. Since waves can travel in the positive and negative z 
directions, a sign convention is introduced for the associated voltages and currents, as shown in Fig. [3T] 
where the voltage vectors of the forward and reflected waves have the same direction and the current 
vectors have opposite directions. 

Using the same time and location dependence as for the electric and magnetic fields in a wave- 
guide, we can write 



/ = \Kkz-CDt) _^\i{-kz-Wt) ^ 

Zq Zq 



(180) 
(181) 



where we have introduced a reflection coefficient P. If we connect a cavity to an impedance at z = 0, 
the expressions above simplify to 



V = VQe-''"'{l+r) 



1= e 
Zo 



:i-r). 



(182) 
(183) 



and the cavity impedance can be expressed in terms of the transmission line impedance Zo and the 
reflection coefficient F: 



^c-y-zo^. 

We can then rearrange the equation for the reflection coefficient and obtain 

P^ Z^-Zo ^ l-jS 
Z^+Zo l+j8' 



(184) 



(185) 




From this equation, we can see that the reflection disappears only for = Zq, the 'matched condition', 
where the waveguide impedance equals the cavity impedance. In the case without beam, this corresponds 
to a coupling factor of j8 = 1 . 

In the context of matching, we therefore have to consider the following points: 

- At the resonance frequency, the power coupler transforms the cavity impedance into the impedance 
of the waveguide. 

- If the cavity is resonating off-resonance or if the coupler is mismatched, power is reflected and 
travels back to the RF source. 

- Since the cavity impedance depends on the Q of the cavity, and since in reality most cavities have 
different Q values, every cavity needs a different matching. 

- Beam loading increases the power needed in the cavity and changes the loaded Q and the cavity 
impedance. Power couplers are usually matched for the case with beam loading. 

- During the start of an RF pulse (before the arrival of the beam), when the cavity is being 'filled' 
with RF power, the cavity is always mismatched, which means we need to make sure that the 
reflected power does not damage the RF source (e.g., by using a circulator between the cavity and 
the RF source). 

The last point is especially important in the case of superconducting cavities, where the dissipated 
power is negligible with respect to the power taken by the beam. In this case one has, basically, full 
reflection of the RF wave at the beginning of the RF pulse before the cavity field increases to its nominal 
level. At that point the beam should enter the cavity, and from then onwards the RF generator is matched 
to the power needs of the RF cavity. After the RF signal is switched off, the cavity voltage decays 



exponentially, as shown in Fig. 32 (more details can be found in 1 18 1). 
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Appendices 

A Cartesian coordinates (x, y, z) 
A.1 Differential elements 



dy 
\dz) 



path element 



dV = dxdydz. 



volume element 



A.2 Differential operators 



V<^> = 



dx 



gradient 



^ dux duy da^ 
dx dy dz ' 



Vxa: 



/ da^ 


day\ 


57 




dux 


_ ^ 


dz 


dx 


day 


dux 


\dx 


dy) 



divergence 



curl 



_ 5^0 

dx^ dy^ dz} 



Laplace 



B Cylindrical coordinates (r, z) 
B.l Transformations 



X = 

y = 



rcoscp, 
rsin<p , 



B.2 Differential elements 



/dr\ 

rd(p 
\dz) 



path element 



(B.2) 



dV = rdrd(pdz- 



volume element 



(B.3) 



B.3 Differential operators 



V<j) 



dr 

1 

r d(p 



gradient 



(B.4) 



Va 



Vxa: 



1 d{rar) 1 da(p da^ 
r dr r d(p dz ' 



/ 1 da^ duq) 



\ 



r dcp dz 

dar da^ 
dz dr 

1 / d{raip) dar\ 
\r \ dr d(p J J 



divergence 



curl 



(B.5) 



(B.6) 



d^d) 1 dd) 1 d^d) d^<b 
ar^ r or r^ a(p'^ dz 



Laplace 



(B.7) 



C Spherical coordinates (r, (p) 
C.1 Transformations 



X = r sin cos 9, 
y = r sin sin 9, 
z = rcost? , 



(C.l) 



with 0<r<oo,o<i><;r,0<<p<2;r. 



C.2 Differential elements 

/ dr \ 
dl = rd7> 

yrsini>d<py 



path element 



dV = sin dr dt? d<p . 



volume element 



C.3 Differential operators 

/ ^ \ 

5r 



1 

1 d(l) 



gradient 



Va = 



1 d{r^ar) 1 (9(a^sint?) 1 da 



+ 



rsinr? 



+ 



Vxa: 



/ 1 /5(a^sint>) 

if 1 dur d{ra(p) 



V 



rsint? ^(p ' 



1 / d{ra^) dttr 



1 



' r dr ' sin d- d-d- d-d- 



+ 



1 



r^sm 



divergence 



curl 



Laplace 



D Useful relationships 



V • (a X b) = b • (V X a) - a • (V X b) , 
V-(Vxa) = 0. 



